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Perturbative analysis of generally nonlocal spatial optical solitons 
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In analogy to a perturbed harmonic oscillator, we calculate the fundamental and some other 
higher order soliton solutions of the nonlocal nonlinear Schodinger equation (NNLSE) in the 2nd 
approximation in the generally nonlocal case. Comparing with numerical simulations we show that 
soliton solutions in the 2nd approximation can describe the generally nonlocal soliton states of the 
NNLSE more exactly than that in the 0th approximation. We show that for the nonlocal case 
of an exponential-decay type nonlocal response the Gaussian-function-like soliton solutions can't 
describe the nonlocal soliton states exactly even in the strongly nonlocal case. The properties of 
such nonlocal solitons are investigated. In the strongly nonlocal limit, the soliton's power and phase 
constant are both in inverse proportion to the 4th power of its beam width for the nonlocal case of 
a Gaussian function type nonlocal response, and are both in inverse proportion to the 3th power of 
its beam width for the nonlocal case of an exponential-decay type nonlocal response. 

PACS numbers: 42.65.Tg , 42.65.Jx , 42.70.Nq , 42.70.Df 
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I. INTRODUCTION 

Since Snyder and Mitchell's pioneering work 1], spa- 
tial solitons propagating in nonlocal nonlinear media 
have been investigated experimentally and theoreti- 
cally in a variety of configurations and material sys- 
tems. It is theoretically indicated that stable spa- 
tial bright (dark) soliton states can be admitted in 
self-focus(self-defocus) weekly nonlocal mediaQ and 
Gaussian-function-like bright soliton states can be admit- 
ted in self-focus strongly nonlocal mediaP,Q . It has been 
shown theoretically that nonlocality drastically modifies 
the interaction of dark solitons by inducing a long-range 
attraction between them, thereby permitting the forma- 
tion of stable dark soliton bound states 4] . The propaga- 
tion properties of light beams in the presence of losses in 
the strongly nonlocal case are different from that in the 
local case[5|. By considering the special case of a loga- 
rithmic type of nonlinearity and a Gaussian function type 
nonlocal response, the dynamics of beams in partially 
nonlocal media 6] and the propagation of incoherent op- 
tical beamsjjj are analytically studied. By using the vari- 
ational principle, the propagation properties of a solitary 
wave in nonlinear nonlocal medium with a power function 
type nonlocal response are studiedQ. The modulational 
instability of plane waves in nonlocal Kerr mediajfj, |l(j 
and the stabilizing effect of nonlocality 01 have been 
studied. The analogy between parametric interaction in 
quadratic media and nonlocal Kerr-type nonlinearities 
can provide a physically intuitive theory for quadratic 
solitons 12J. Some properties of the strongly nonlocal 
solitons (SNSs) and their interaction are greatly differ- 
ent from that in the local case, e. g. two coherent SNSs 
with 7T phase difference attract rather than repel each 
other 0, the phase shift of the SNS can be very large 
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comparing with the local soliton with the same beam 
width[j|, and the phase shift of a probe beam can be 
modulated by a pump beam in the strongly nonlocal 
case|l3j. Employing a Gaussian ansatz and using a vari- 
ational approach, the evolution of a Gaussian beam in 
the sub-strongly nonlocal case is studied 14]. Recently 
it is experimentally shows that solitons in the nematic 
liquid crystal(NLC) are SNSs [H HI]. The team of As- 
santo has developed a general theory of spatial solitons 
in the NLC that exhibiting a nonlinearity with an arbi- 
trary degree of an effective nonlocality and established 
an important link between the SNS and the parametric 
soliton^, ITfl IItI Hsj . They also experimentally inves- 
tigated the role of the nonlocality in transverse modu- 
lational instability(MI) in the NLC [lit and observed 
the optical multisoliton generation following the onset of 
spatial Ml[l9j. The interaction of SNSs has been ex- 
perimentally demonstrated [2fl l2l| . and the possibility of 
all-optical switching and logic gating with SNSs in the 
NLC has been discussed [22^ . 

However the theoretical studies on the spatial nonlo- 
cal soliton are mostly focused on the strongly nonlocal 
casep], 0, IE 0, EI \JM an d the weekly nonlocal case0- 
There is a lack of studies on the moderate nonlocal case. 
On the other hand, even though a convenient method 
has been introduced in references [3, la, Hia, to study 
the propagation of light beams in the strongly nonlocal 
case or even in the sub-strongly nonlocal case, to em- 
ploy this method efficiently the nonlocal response func- 
tion must be twice differentiable at its center. As will 
be shown this method can't deal with the nonlocal case 
of an exponential-decay type nonlocal response function 
that is not differentiable at its center. In this paper, 
in analogy to a perturbed harmonic oscillator, we calcu- 
late the fundamental and some other higher order soli- 
ton solutions of the NNLSE in the 2nd approximation in 
the generally nonlocal case. Our method presented here 
can deal with the nonlocal case of an exponential-decay 
type nonlocal response function. Numerical simulations 
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conform that the soliton solution in the 2nd approxima- 
tion can describe the generally nonlocal soliton states of 
NNLSE more exactly than that in the Oth approximation. 
It is shown that for the nonlocal case of an exponential- 
decay type nonlocal response the Gaussian-function-like 
soliton solutions can't describe the fundamental soliton 
states of the NNLSE exactly even in the strongly nonlo- 
cal case, that is greatly different from the case of a Gaus- 
sian function type nonlocal response. The properties of 
such nonlocal solitons are investigated. The functional 
dependence of such nonlocal soliton's power and phase 
constant on its beam width is greatly different from that 
of the local soliton. Further more this functional de- 
pendence for the nonlocal case of a Gaussian function 
type nonlocal response greatly differs from that of an 
exponential-decay type nonlocal response. In particu- 
lar in the strongly nonlocal limit, the nonlocal soliton's 
power and phase constant are both in inverse proportion 
to the 4th power of its beam width for the nonlocal case 
of a Gaussian function type nonlocal response, and are 
both in inverse proportion to the 3th power of its beam 
width for the nonlocal case of an exponential-decay type 
nonlocal response. 

II. THE FUNDAMENTAL GENERALLY 
NONLOCAL SOLITON SOLUTION IN THE 2ND 
APPROXIMATION 

Let's consider the (1+1)-D dimensionless nonlocal non- 
linear Schodinger equation(NNLSE) 0, H B H ES [13 

8v 1 8 2 v r +OQ 

l Tz + 2d^ + u J „ fl(*-OK^)l 2 ^ = o, (i) 

where u(x, z) is the complex amplitude envelop of the 
light beam, x and z are transverse and longitude co- 
ordinates respectively, R(x) > is the real symmetric 
nonlocal response function, and 

/+oo 
R(x-£)\u(tz)\ 2 dt; (2) 
-oo 

is the light-induced perturbed refractive index. 

As indicated in reference 3], if R(x) is twice differ- 
cntiable at x — and the 2nd derivative R"(0) < 0, 
and if the characteristic nonlocal length is one order of 
the magnitude larger than the beam width of the soli- 
ton, the NNLSE Q can be simplified to the following 
strongly nonlocal model (SNM) 

i Pz + 1 2^ +U J +00[Ro + % X - Z)|2rf? = °' 

(3) 

where R = R(0) and R'^ = R"(0). For example, for 
the Gaussian function type nonlocal response function 
R(x) = 1/(w \/tt) exp (— x 2 /wq) , when the character- 
istic nonlocal length wq is one order of the magnitude 



larger than the beam width, the SNM J2J can describe the 
NNLSE very well 3J. However, as will be shown, when 
the characteristic nonlocal length and the beam width are 
in the same order of the magnitude, the SNM can't 
describe the NNLSE JTJ very well. The SNM © can't 
deal with the generally nonlocal case. Further more, for 
the exponential-decay type nonlocal response function 
R{x) — l/(2?«o) exp(— l^l/uio) which is not differentiable 
at x = 0, we can't get the parameter Rq of the SNM J3Jl. 
So the SNM can't deal with this nonlocal case of such 
an exponential-decay type nonlocal response. 

The SNM © allows a Gaussian-function-like bright 
soliton solution 

/ 1 \V4 x 2 3 

u (x, z) = A (— cxp [___ J (__ i?0 A 2 )z] (4) 
where 

^ = -KA\ (5) 

and v is the beam width of uq(x, z). The power P and 
the phase constant 7 of uo(x, z) are given by 

/+00 
\u(x,t)\ 2 dx = A 2 : (6) 
-OO 

7 = RoA* - (7) 

respectively. 

In this paper, we define the degree of nonlocality by the 
ratio of the characteristic nonlocal length to the beam 
width of the light beam and use the phrase "generally 
nonlocal case" to refer to the nonlocal case where the 
degree of nonlocality is larger than one and less than 
ten. For the Gaussian function type nonlocal response 
function and the soliton solution (0} , the degree of non- 
locality is wq/v. The larger of wo/v, the stronger of the 
nonlocality. In fact for a given type of nonlocal response, 
soliton solutions with the same degree of nonlocality can 
be described in the same way. That can be clarified by 
taking transformations pijj 

x — — z — — u — ku R — kR . (8) 

K K 

Under these transformations, the form of NNLSE (JTJ 
keeps invariant and the degree of nonlocality keeps in- 
variant too. If we set n equal to the characteristic nonlo- 
cal length of R(x), the characteristic nonlocal length of 
R(x) will be scaled to unity and the degree of nonlocal- 
ity will be determined only by the beam width of u(x, z). 
In this case the less of the beam width of u(x,z), the 
stronger of the nonlocality. On the other hand we may 
also set K equal to the beam width of u(x, z). If we do 
this, the degree of nonlocality will be determined only by 
the characteristic nonlocal length of R(x). The larger of 
the characteristic nonlocal length of R(x), the stronger 
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of the nonlocality. In this paper, the characteristic non- 
local length of R(x) and the beam width of u(x, z) are 
not scaled to unity. 

For the soliton state u(x,z), we have \u(— x,z)\ 2 — 
\u(x, z)\ 2 and \u(x, z)\ = |tt(x,0)|. So for the soliton state 
u(x,z), by defining 



V(x) = - 



+ 00 



R(x-OH^z)\ 2 d^, 



the NNLSE JTJ reduces to 
.du 1 d 2 u 



dz 



2dx 2 



V{x)u = 0. 



(9) 



(10) 



Taking the Taylor's expansion of V(x) at x — 0, we ob- 
tain 



where 



V{x) =V + ^-rx 2 + ax 4 + Px 6 + 
2/i 4 



V = V(0), 
^ = V^>(0), 

a=£V< 4 >(0), 



(11) 



(12a) 
(12b) 

(12c) 

(12d) 



As will be shown, in the generally nonlocal case and the 
strongly nonlocal case the parameter u can be viewed 
as the beam width of the soliton, and when x < /i, the 
terms ax 4 and f3x e are one and two order of the magni- 
tude smaller than the term x 2 /(2u 4 ) respectively. That 
indicates the effects of ax 4 and (3x & on the soliton are 
considerably small comparing with the effect of x 2 /(2u 4 ) 
in the generally nonlocal case. Further more in the gen- 
erally nonlocal case, the effects of the a; 8 power term 
and the other higher power terms of the Taylor's series 
of V(x) on the soliton are far smaller than the effects 
of these three lower power terms. For convenience sake 
we will neglect such higher power terms in the following 
discussions and simply adopt 
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V(x) = V + -^x z + ax 4 + (3x°. 



(13) 



However, as the degree of nonlocality decreases the effects 
of ax , (3x 6 and other higher power terms become larger 
and larger, and when the characteristic nonlocal length 
is comparable with or less than the beam width of the 
soliton, the x 8 power term and other higher power terms 
are no longer negligible. For such cases we must take the 
higher power terms of the Taylor's series of V(x) into 
account. 

In the generally nonlocal case, substitution of Eq. 113(1 
into Eq. (TDJ yields 



du 
oz 



1 d 2 

20X 2 



To 



1 



/3x e 



u. (14) 



Taking a transformation 



u(x, z) = ip n {x) exp[-i(e„ + Vo)z] 



we arrive at 
2dx 2 



2/i 4 



Px t 



(15) 



(16) 



where the index n = 0, 1, 2, • • • is the order of the soliton 
solution, in particular n — corresponding to the funda- 
mental soliton solution and n = 1 corresponding to the 
second order soliton solution and so on. Even though 
Eq. ((16(1 takes the form of the stationary Schrodinger 
equation, the parameters [i, a, (3 are depended on the soli- 
ton solution ip n (x). 

If a = and (3 = 0, equation ((16(1 reduces to the 
well-known stationary Schrodinger equation for a har- 
monic oscillator. Since in the generally nonlocal case the 
effects of the terms ax 4 and (3x & on the soliton are far 
smaller than that of the term x 4 /(2fj, 4 ), we view the terms 
ax 4 and f3x 6 as perturbations in the process of solving 
Eq. 1 (16(1 . Following the perturbation method presented 
in any a text book about quantum mechanics (for exam- 
ple, seeing |24|), we get for the fundamental soliton solu- 
tion in the 2nd approximation 



MA a, /?, x) « A{^f' 4 exp(-^) 
x[l + a( 9 4 
+a 2 ( 



16 

j , 1247m 1 - 
512 



3^ x 2 _ Ml ^4 



4 

141m 1 
64 



^ 4 ) 



53/i 
64 



13m ~6 , lS~S\ 
48 ~ 32 I 



a/55M_ _ lajjr 2 _ an' x 4 _ M_y 
' y 32 8 8 6 A 



15m 



■ 



(17) 



and 



so 



2(i 



1 3^ 4 a 21y}°a 2 Wu 6 P , . 
- + ^ + ^- (18) 
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In Eqs. 1(171) and (|18[). if we neglect the a, a 2 , (3 terms 
or neglect the a 2 , (3 terms only, we will get for the fun- 
damental soliton solution in the 0th approximation or in 
the 1st approximation respectively. 

Substituting Eq. (TJJ into Eq. © , we have 



V(A,a,0,x) 



+00 



R{x-£)il$(A,a,0,t)dt. (19) 



Keeping in mind Eqs. 1(12(1 . we obtain 

= y( 2 )(A,a,/3,0), 
= iy( 4 )(A«,/3,0), 
/?= |^ 6 )(A,a,/?,0). 



A* 

a - 



(20a) 
(20b) 
(20c) 



For a fixed value of parameter /Lt, the parameters A, a, [3 
can be found by solving Eqs. ((20(1 . In the section Ap- 
pendix (|A")l . we present a fixed-point method to numer- 
ically calculate these parameters based on Eqs. ((20(1 . 
Hereinabove we have formally presented the main for- 
mulas to calculate the perturbed fundamental generally 
nonlocal soliton solution in the 2nd approximation. 



4 



A. the nonlocal case of a Gaussian function type 
nonlocal response 

As an example, let us consider the nonlocal case of a 
Gaussian function type nonlocal response pi Irl PI Wi Hlj 



R(x) 



exp 



(21) 



For the SNM © and the soliton solution (QJ, we can 
find the fundamental soliton solution for such a Gaussian 
function type nonlocal response in the strongly nonlocal 



u (x, z) = (- 



2^ 



4 e 



(22) 



This soliton solution can describe the soliton state of the 
NNLSE JIJ exactly in the strongly nonlocal case when the 
degree of nonlocality w^/v > 10, but can't describe the 
soliton state in the generally nonlocal case when wq/v ~ 
2. 

In the generally nonlocal case, the fundamental soli- 
ton solution in the 2nd approximation is described by 
ipo(A,a, 0,x) in Eq. ((T7|l . As shown in Fig. QJ when 
w a = 2, fi = 1 and A = 3.22, a = -0.0487, = 
0.00317 numerically calculated by the fixed-point method 
presented in the section Appendix (^J, the difference 
between the fundamental soliton solution in the 2nd 
approximation ipo(A,a,0,x) and that in the 0th ap- 
proximation tpo(A, 0, 0, x) is comparatively small. As a 
Gaussian function, the power and the beam width of 
ij)o(A,Q, 0,x) are given by A 2 and fi respectively. There- 
fore the power and the beam width of ipo(A, a,0, x) are 
approximatively given by A 2 and [i respectively too. So 
in the generally nonlocal case we can approximately de- 
termine the degree of nonlocality by wo/[x, and approxi- 
mately obtain 



A 2 



exp 



exp 



^(A,0,0,£R 



and 



V 



a •■ 

0'^ 



^(i+ w 2 /„ 2 )^ 2 

2ju ' 
_ (1+wl/v 2 ) 

1 

4^ (1+w 2 /^ 2 ) ' 

1 

12^(1+^/^2)2 ■ 



(23) 

(24a) 

(24b) 
(24c) 
(24d) 



Using Eqs. I|24(l for wq — 2 and /j, = 1, we can find 
A sw 3.14, a ps — i and s» that are very close to 
the numerically calculated values A = 3.22, a = —0.0487 
and = 0.00317, and we can find | x y^) i ~ liS^I < 

?\ < 0.007 for x < n that are 



0- 1 and l xV(V) 



I 150m 4 



In the strongly nonlocal limit the degree of nonlocality 
w o/m ^ 1) we have 



A 2 

Va 

a rs 
/3« 



2m 4 ' 
2m 3 "' 



(25a) 
(25b) 
(25c) 
(25d) 



As the degree of nonlocality wq/ijl approaches in- 
finity, the parameters a, both approach zero, and 
ipo(A, a, 0, x) approaches ipo(A, 0, 0, x). In such a case a 
Gaussian-function-like strongly nonlocal soliton solution 
is obtained. 

Using the NNLSE 0) as the evolution equation 
and using the numerical simulation method we investi- 
gate the propagation of light beams in nonlocal media 
with a Gaussian function type nonlocal response. The 
numerical simulation method is the split-step Fourier 
Method(SSFM)|2j|, the step-size Az = 0.01, transver- 
sal sampling range —10 < x < 10 and the sampling 
interval Ax = 0.1. With different input amplitude en- 
velops (the initial data of numerical simulations) that 
are described by uq(x,0) in Eq. {23)1. ipo(A,0,Q,x) and 
ij)o(A,a,0,x) respectively, we show the propagations of 
these light beams in Fig. (0. It is indicated that in the 
generally nonlocal case when the degree of nonlocality 
Wq/h = 2, (A, a, 0, x) can describe the soliton state 
of the NNLSE more exactly than i/; (A, 0, 0, x) and 
uq{x,Q) in Eq. (|22|) . The soliton solution in the 2nd ap- 
proximation iI)q{A, a, 0, x) also can describe the soliton 
state of the NNLSE exactly when Wo/fi = 1, that is 
shown in Fig. (|3J). However when wq/^i = 0.5, as indi- 
cated in Fig. ipo(A, a, 0, x) can't describe the soliton 
state of the NNLSE exactly. In such a case, we must 
take the higher power terms of the Taylor's series of V(x) 
into account and calculate the higher order approxima- 
tion. To show how exactly i^o{A, a, 0, x) describe the 
fundamental soliton state, we define 



9{z) 



\u(x,z)e~ i 't > ( z )-u(x,0)\ 2 dx 
J+™\u(x,0)\ 2 dx ' 

E _ /' e(z)dz 



(26a) 
(26b) 



where e 1 ^^ is the phase factor of u(x, z) and for the 
fundamental soliton e 1 ^ 7 -) = ^j^ri ■ For a fixed value of 

I, the less of 8, the more exactly ipo(A, a, 0, x) describe 
the fundamental soliton state. As shown in table |IJ), 
tpo(A, a, 0, x) can describe the fundamental soltion states 
exactly when Wo/fi > I. 

Now let us consider the properties of ipo(A,a,0,x), 
As have been shown, the beam width of ipo(A,a, 0,x) 
is approximatively given by fi, and its power and phase 
constant are approximatively given by 



consistent with the perturbation postulate. 



V¥(l + w 2 /^/ 2 
2fi 



(27) 



FIG. 1: The comparison between \ipo(A, a, P,x)\ 2 (dashing 
line) and \tpo{A, 0, 0, x)\ 2 (solid line). Here wo = 2, fi — 1, A — 
3.217, a = -0.0487, (3 = 0.00317 and the degree of nonlocality 
Wo/n = 2. 



TABLE I: Using the numerical simulation method we cal- 
culate 9 in Eqs. 1261 for the nonlocal case of the Gaussian 
function type nonlocal response and for the nonlocal case of 
the exponential-decay type nonlocal response. 



e 




0.020 c 0.0033 d 
0.044 c 0.013 d 


0.0027 c 0.0027' 0.0026 9 
0.013 E 0.013 / 0.013 9 


0.0026 h 
0.013' 1 


0.0026* 
0.013 s 


4>o b 
V2 6 


0.017 c 0.0095 d 
0.072 c 0.044 d 
0.36 c 0.20 d 


0.0078 c 0.0072' 0.0071 9 
0.037 e 0.034 / 0.032 s 
0.17 e 0.15 / 0.14 9 


0.0069" 
0.031 h 
0.14' 1 


0.0068 8 
0.031* 
0.13 s 



"the nonlocal case of the Gaussian function type nonlocal response 

6 the nonlocal case of the exponential-decay type nonlocal 

c thc degree of nonlocality equal to 1 

d the degree of nonlocality equal to 2 

e the degree of nonlocality equal to 3 

'the degree of nonlocality equal to 4 

'the degree of nonlocality equal to 5 

'the degree of nonlocality equal to 6 

'the degree of nonlocality equal to 7 



7 = -Va - £0 



1 



8(l+™g/M 2 ) " + " 64(1+™ Hp?Y 



(28) 



respectively. In the strongly nonlocal limit the degree of 
nonlocality Wq/ /j> 1, we have 



P 



2 M 4 : 



2[x 



4 • 



(29) 



(30) 



That means for a given value of characteristic nonlocal 
length in the strongly nonlocal case the power and the 
phase constant of the nonlocal soliton are both in inverse 
proportion to the 4th power of its beam width. The 
dependence of the power P and the phase constant 7 
on the beam width /1 are shown in Fig. JSJ for a given 
value of characteristic nonlocal length. It is indicated 
that Eq. l(2T|) and Eq. (|3Tjj) can describe these dependence 
exactly in the generally nonlocal case. 

To make a comparison with the local soliton, let 
us consider the following local nonlinear Schrodinger 






FIG. 2: The propagations of light beams in nonlocal media 
with a Gaussian function type nonlocal response with differ- 
ent input intensity profiles that are described by (a) \uo(x, 0)| 2 
in Eq. Q3, (b) \ipo(A, 0, 0, x)\ 2 and (c) \ipo{A, a, fi, x)\ 2 re- 
spectively. Here too = 2,/^ = 1, v = I, A = 3.217, a = 
-0.0487, P = 0.00317 and the degree of nonlocality w /m = 2 - 



equation(NLSE)|2J 

.du 

az 



ld 2 u 
2dx^ 



+ \u\ 2 u = 0. 



(31) 



When the characteristic nonlocal length wq approaches 
zero, the Gaussian function type nonlocal response 
function R(x) approaches the 6(x) function, and the 
NNLSE (JTJ approaches the NLSE J^D- The fundamental 
soliton of the NLSE (JSTJ) is given bv|26| 



u(x, z) — — sech(— ) exp ( i— ^ 
rj r] \ 2r\ jL 



(32) 



where 77 can be viewed as the beam width of the local 
soliton. The power and the phase constant of such a 
local soliton are given by 



P 



+00 2 
\u(x, t)\ 2 dx = —. 
■00 V 



(33) 



6 




FIG. 3: (a) The comparison between R(x) and 
\ipo{A, a, [5, x)\ 2 . Dashing line: R(x)/R(0); solid line: 
\ipoiA,a,P,x)f/\ipo(A,a,P,0)\ 2 ; (b)The propagation of 
the light beam with an input intensity profile described 
by \ip (A,a,/3,x)\ 2 . Here w = I, (J, = I, A = 1.777, a = 
—0.113, (3 = 0.0177 and the degree of nonlocality wo/n = 1. 



FIG. 4: (a) The comparison between R(x) and 
\ipo{A,a, f3,x)\ 2 . Dashing line: R(x)/R(0); solid line: 
\ipo{A,a,P,x)\ 2 /\i[)o(A,a,/3,0)\ 2 ; (b)The propagation of the 
light beam with an input intensity profile described by 
\ipo(A,a,p,x)\ 2 . Here w = 0.5, n = 1,A = 1.403, a = 
— 0.178, p — 0.0476 and the degree of nonlocality Wo/fi = 0.5 



case, equation (|33|l leads 



1 



(34) 



respectively. We can find the power and the phase con- 
stant of the local soliton are in inverse proportion to the 
1st and the 2nd power of its beam width respectively. 
The functional dependence of the power and the phase 
constant of the nonlocal soliton on its beam width greatly 
differs from that of the local soltion. 



B. the nonlocal case of an exponential-decay type 
nonlocal response 

As another example, we investigate the nonlocal 
case where the light-induc ed p erturbed refractive index 
n(x, z) is governed by|j, Ell 



n{x, z) 



wl 



d 2 n(x, z) 
dx 2 



n{x,z) = - — 
2w 



-Hoc. 



exp 



l*-g| 
w 



|«&*)| 2 de, (36) 



and we get the exponential-decay type nonlocal 
response HBSIISini 



R(x) = — expf-M 
2wo V wo 



(37) 



Since the exponential-decay type nonlocal response func- 
tion R(x) is not differentiable at x = 0, the SNM © 
can't deal with this nonlocal case. So we have to 
use ipo (A, a, /3, x) to describe the soliton state of the 
NNLSE III). 

For this exponential-decay type nonlocal response and 
the fundamental soliton state, V(x) can be approxi- 
mately given by 



\u(x,z)\ 2 =0. (35) V(x)fn-J_ 



+oo i 
oo 2wo 



exp 



It is found that several nonlocal media, for example the 
nematic liquid crystal [THl ITsj . their light-induced per- 
turbed refractive index can be described by Eq. (|35|l . If 
the size of the nonlocal media is much larger than the 
beam width of the soliton and the characteristic nonlo- 
cal length, the effect of the boundary condition on the 
soliton can be negligible and we can simply assume the 
size of the nonlocal media is infinity large. For such a 



4wq ^ ° ° y 2wq fi 



i] 



- e ™° l er f 



where 



erf(x) 



1]}, 



(38) 



(39) 
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FIG. 5: The dependence of the power P and the phase con- 
stant 7 on the beam width fi. Dashing lines are described 
by Eq. I|27jl and Eq. i'AOt respectively; Thin solid line is de- 
scribed by Eq. 1291 : Thick solid lines are directly calculated 
with parameters A, a, (3 numerically found by the fixed-point 
method presented in the section Appendix 10. Here wo = 6. 



FIG. 6: The propagations of light beams with different input 
intensity profiles that are described by (a) \ipo{A, 0, 0, x)\ 2 and 
(b) \%po(A, a, f3, x)\ 2 respectively. Here wo = 100, /i = 1,A = 
138.159, a = -0.0767, (3 = 0.0150 and the degree of nonlocal- 
ity wo/fJ. = 100. 



Combining Eqs. (|T2l . we get 



y» 



Vn 



48 ir 



1440tu,( 



(40a) 

(40b) 
(40c) 
.(40d) 



In the strongly nonlocal limit the degree of nonlocality 
wo/^^ 1, we obtain 



A 2 



V 



\/™0 

^ ' 

2^ ■ 
1 

12/j 6 ' 
1 

60p 8 ' 



(41a) 
(41b) 
(41c) 
(41d) 



It is worth to note that in the strongly nonlocal case the 
parameters a and f3 are free from the characteristic non- 
local length wq. Even when the characteristic nonlocal 
length wq approaches infinity, the parameters a,/3 still 
rest on finite values and don't approach zero, and there- 
fore tpo(A,a, f3,x) does't approach il>o{A,Q,0,x). That 
greatly differs from the nonlocal case of a Gaussian func- 
tion type nonlocal response. As a result the Gaussian- 
function- like soliton solution ipo (A, 0, 0, a;) can't describe 
the soliton state of the NNLSE (JTJ exactly even in the 



strongly nonlocal case, that is shown in Fig. JBJ. As 
shown in Fig. J7|), ipo(A, a, (3, x) also can describe the soli- 
ton state of the NNLSE exactly when u>o//j = 1. Even 
when wq/h = 0.5, as indicated in Fig. ©, i^^A^a, (3,x) 
can describe the soliton state of the NNLSE in high 
quality. As indicated by the values of 8 in table {U, 
ipo(A, a, /3, x) can describe the fundamental soliton states 
of NNLSE Q exactly in the generally nonlocal case. 

In the strongly nonlocal case the soliton's power and 
phase constant are approximately given by 



Pa A 2 



-Vn 



/TTWn 



2^ 3 



(42) 



(43) 



respectively. For a given value of the characteristic non- 
local length, the soliton's power and phase constant are 
both in inverse proportion to the 3th power of its beam 
width in the strongly nonlocal case that differs from the 
nonlocal case of a Gaussian function type nonlocal re- 
sponse where the soliton's power and phase constant are 
both in inverse proportion to the 4th power of its beam 
width in the strongly nonlocal case. The dependence of 
the soliton's power P and phase constant 7 on its beam 
width /i are shown in Fig. @ for a given value of char- 
acteristic nonlocal length. It is indicated that Eq. I|42|) 
and Eq. I|43|l can describe these dependence very well in 
the strongly nonlocal case. 
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FIG. 7: (a) The comparison between R(x) and 
\ipo{A,a, f3,x)\ 2 . Dashing line: R(x)/R(0); solid line: 
\ipo(A,a,f3,x)\ 2 /\ipo(A,a,f3,0)\ 2 ; (b)The propagation of 
the light beam with an input intensity profile described 
by \tpo{A, a, /3,x)\ 2 . Here w = = I, A = 2.206, a = 
—0.126, /? = 0.0280 and the degree of nonlocality wo/fJ,— 1. 



FIG. 8: (a) The comparison between R(x) and 
\ipo{A,a, f3,x)\ 2 . Dashing line: R(x)/R(0); solid line: 
\ipo(A,a, /3,x)\ 2 /\ij)a{A,a, /3,0)| 2 ; (b) The propagation of 
the light beam with an input intensity profile described 
by \ip (A,a,/3,x)\ 2 . Here w = 0.5, n = 1,A = 1.610, a = 
—0.158, j3 — 0.0408 and the degree of nonlocality wo/m = 0.5. 



III. THE HIGHER ORDER GENERALLY 
NONLOCAL SOLITON SOLUTIONS IN THE 2ND 
APPROXIMATION 

A. the nonlocal case of the Gaussian function type 
nonlocal response 

The second order soliton solution for the SNM J3J is 
given by 



the 2nd approximation is given by 



^ [A, a, /?, x) » A(^) 1 / 4 exp(-|l) 



[l + a(^f-^fx 2 -^ 



21L. 



+Q 2 (- 



+/?( 



8375^ 
512 

385p 8 



32 



215m 1u r 2 
64 ^ 



3V1 7 ,2 
8 



73^ 
64 



19f ° ~6 , Hl^8l 
48 ~ 32 I 



(47) 



Ui(a;,z) = A ( — ^ 

1 -KV £ 



where 



1/4 



V2a 



-e 2» 



- i (ifV-#o J 4 2 )z 



. (44) 



and 



si 



3 



15a/i 4 165aV° 



105/?// 



(48) 



" /l 2 



(45) 



The power and the beam width of u\(x,z) are given by 
A 2 and \fiv respectively. For the Gaussian function type 
nonlocal response, the second order soliton solution for 
the SNM Q is given by 



ui(x,z) = ( 



2v 4 



e 2,/ 



(46) 

This soliton solution can describe the second order soliton 
state of the NNLSE exactly in the strongly nonlocal 
case when wq/{^/%v) > 10 but can't describe it exactly 
in the generally nonlocal case when wq/{\/Zv) ~ 2. 
The second order generally nonlocal soliton solution in 



For the Gaussian function type nonlocal response func- 
tion l |2ip . as shown in Fig. I)l()[l. the difference between 
l^i (A /3, a;)| 2 and \ipi(A, 0, 0, x)| 2 is small in the gen- 
erally nonlocal case. As an Hermite-Gaussian function, 
the power and the beam width of ipi(A, 0, 0, a;) are given 
by A 2 and \/3m respectively. So the power and the beam 
width of i/ji(A, a, (3, x) are also approximatively given by 
A 2 and V3/i respectively. We can approximatively de- 
termine the degree of nonlocality by wo/(V3fi) an( i a P~ 
proximatively abtain 



V(x)*>-£ 



exp 



5 — 



^(A,0,0,0« 



•7j2+^ 2a; 2 / A1 2 + tt .g/ A1 2 + t 
(1+^/M 2 ) 2 



4 / 4 



(49) 
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1 1.5 2 2.5 3 3.5 4 
beam width 

FIG. 9: The dependence of the soliton's power P and phase 
constant 7 on its beam width /j,. Dashing lines are described 
by Eq. 1 12t and Eq. IJ430 respectively; Solid lines are directly 
calculated with parameters A, a, (3 numerically found by the 
fixed-point method presented in the section Appendix lIXl . 
Here wq — 6. 




-2 - 1 1 2 



FIG. 10: The comparison between \ipi(A,a, f3,x)\ 2 (dashing 
line) and l^i (^4, 0, 0, x)\ 2 (solid line). Here wo = 1.5, (1 = 
l/\/3, A = 7.606, a = -0.450, (3 = 0.00530 and the degree of 
nonlocality u;o/(\/3/i) = 1.5. 



a 




-4-2 2 4 



and 



.4 2 



2\5/2 



2) ' 

_ ^(l+u-g/p 2 ) 
~ 2^(^/ M 2-2)' 

K/^-4) 
4 M 6 (l+^/M 2 )K/M 2 -2)' 



K/M 2 -6) 

12^(1 + ^/^2)2(^/^2 



(50a) 
(50b) 
(50c) 
(50d) 



In the strongly nonlocal limit the degree of nonlocality 
Wo/(V3fi) » 1, we have 



A 2 

V 

a ps 



J 2^i 4 ' 

2/i 4 ' 
1 

4/j 4 io(; ' 
1 

12p 4 M) 4 ' 



(51a) 
(51b) 
(51c) 
(51d) 



As the degree of nonlocality u>o/(v3/-t) approaches in- 
finity, the parameters a and /3 approach zero, and 
■01 (A ce, (3, x) approaches ipi(A, 0, 0, x). Therefore in the 
strongly nonlocal case an Hcrmitc-Gaussian-function-likc 
second order soliton solution is obtained, and the power 
and the phase constant of ipi(A, a, (3, x) are both in in- 
verse proportion to the 4th power of its beam width. As 
indicated in Fig. (|11|) and Fig. I|12|) . the second order soli- 
ton solution in the 2nd perturbation ipx{A, a, (3, x) can 
describe the second order soliton state of the NNLSE Q 
exactly when u>o/(\/3/i) = 1.5 and describe it in high 
qulity when u>o/(\/3m) — 1- As shown by the values of 




FIG. 11: (a) The comparison between R(x) and 
\ifji(A, a, /?, x)\ 2 . Dashing line: R(x)/R(0); solid line: 
\^i(A,a, f3,x)\ 2 /\il>i(A,a, /3,0.65)| 2 ; (b)The propagation of 
the light beam with an input intensity profile described 
by |i/)i(A,a,/3,a;)| 2 . Here w = 1.5,/i = 1/V3,A = 
7.606, a = -0.450, (3 = 0.00530 and the degree of nonlocality 
u>o/(V3iJ.) = 1.5. 



6 in table $,ipi(A, oc, (3, x) can exactly describe the sec- 
ond order nonlocal soltion state in the generally nonlocal 
cases. 

Finally since the all eigenfunctions of the harmonic os- 
cillator can be found systematically [24|, it is possible that 
in analogy to a perturbed harmonic oscillator we can also 
approximately calculate the third order soliton solution 
or the fourth order soliton solution and so on in the gener- 
ally nonlocal case for the Gaussian function type nonlocal 
response. 
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FIG. 12: (a) The comparison between R(x) and 
\i)ji(A, a, (3, x)\ 2 . Dashing line: R(x)/R(0); solid line: 
\^i(A,a, f3,x)\ 2 /\-il)i(A,a, f3,0.65)\ 2 ; (b)The propagation of 
the light beam with an input intensity profile described by 
\ipi(A,a,f3,x)\ 2 . Here w = 1,/i = l/y/3,A = 6.133, a = 
-0.492, j3 = 0.0145 and the degree of nonlocality w /(V3n) = 
1. 



B. the nonlocal case of the exponential-decay type 
nonlocal response 

As have been indicated, to [A, a, 0, x) and 
ipi(A, a, (3, x) in the nonlocal case of the Gaussian func- 
tion type nonlocal response and to ipo(A, a, (3, x) in the 
nonlocal case of the exponential-decay type nonlocal re- 
sponse, we have V^(A, a, (3, 0) > for generally nonlo- 
cal cases. If we define 



V(x) = V + -^-rx 2 + ax A + (3x 6 , 



(52) 



we will get V(A, a, f3, x) w V(x). But as shown in 
Figs. l(T5k1 and lJT4V l. to (A, a, (3, x) in the nonlocal 
case of the exponential-decay type nonlocal response, we 
have (A, a, (3, 0) < 0. In such a case we can't de- 
fine l//i 4 = V' 2 ) (A, a, (3, 0) and can't define the param- 
eters n,a,(3 as those in Eqs. (|20|l . However as shown in 
Figs. l[T5k) and ifHkl. we still can find suitable values 
of A, a, (3 for a fixed value of \x to make V(A, a, (3, x) « 
V(s). These suitable values of A,a,f3 can be calculated 
by solving the following coupling equations 



V{A,a,j3,x ) = V + 



1 „.2 



+ axQ + /3xg , (53a) 



V'(A, a, f3, x Q ) = ±x + ^axl + 6/3x1 (53b) 



V"(A,a,0,x o ) 



+ I2ax% + 30/3x$, (53c) 



where x ^ and Vo = V(vl, a,(3,0). In the sec- 
tion Appendix |Bj) we present a fixed-point method 



to calculate these parameters A, a, f3 with Eqs. I|53|l . 
In Figs. |^23),|J2t>) an d (fHjb 'l we show the propaga- 
tion of lights with input intensity profiles described by 
|^i (^4, a, f3, x)\ 2 . Even when the wo/(\/3[i) = 0.5, there 
still exists a second order nonlocal soliton. As shown by 
the values of 9 in table {U, ipi (A, a, j3, x) can describe 
the generally nonlocal soliton state in high quality. Since 
the difference between tpi(A, a, {3, x) and il>i(A, 0, 0, x) is 
small, we can approximately get 



^(^+2)e^{e- 



wo 



+e^[er/(^- + £) 



|^(A,0,0,C)| 2 ^ 

1]} 



2-y/Tru 



By defining 

tffc) = V(x)/A 2 
and combining with Eqs. I|53|) . we obtain 



(54) 



(55) 



.4: 



4x 2 / M 4 



24[[/(zo) - U(0)] - 9U'{x )x + U"{x Q )x 2 ' 

(56) 

For example, when wq — 10, fx = l/y3,Xo = 2, from 
Eqs. (EH, and we get A w 47.323 that is close 
to the numerically calculated value A — 48.257. 

While V(A,a, (3,x) w V"(a:), as shown in Figs, ifTSk ) 
and |14k) there still exists difference H(x) = 
V(A, a, p, x) — V(x). To achieve higher accuracy we 
should take H{x) into account and set V{x) — Vo + 
x 2 /(2fi 2 )+ax 4 +j3x 6 + H(x). Viewing H( x ) as perturba- 
tion we will obtain another higher accurate second order 
soliton solution. However the form of H (x) is rather com- 
plex and we will leave it in future further work and don't 
intent to deal with the effect of H(x) in this paper. 

Now let us consider the third order nonlocal soliton. 
The third order generally nonlocal soliton solution in the 
2nd approximation is given by 



24587^ 



45m° 
16 

10 



123^ 



v 7T/i z 

^2 



^)+a 2 C-l^ 



256 **• 

655m 8 
32 



8 

4lM 8 . 
64 ' 

1405m 6 



13m 
4 



X 



193^ 6 
96 ^ 

_ 125m 4 ^4 



97m „.8 , t£_~W\ 
96 ^ 16 i 



25m 
12 



X 



±* 8 )J57) 



16 8 12 3" 

As shown in Fig. (|To|> and table ^,if)\{A,a, (3,x) can 
describe the third order generally nonlocal soliton only 
qualitatively. To obtain a higher accurate third order 
soliton solution we should take all perturbation into ac- 
count or develop another new method. 



IV. CONCLUSION 

In analogy to a perturbed harmonic oscillator, we cal- 
culate the fundamental and some other higher order soli- 
ton solutions in the 2nd approximation in the generally 
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FIG. 13: (a) The comparison between V(A, a, j3, x) (solid 
line) and V(x) (dashing line); (b)The propagation of the 
light beam with an input intensity profile described by 
\ipi{A,a,P,x)\ 2 . Here w = 2,/j, = 1/^3, A = 12.162, a = 
-0.442, p = 0.0179 and the degree of nonlocality too/(%/3ju) = 
2. 
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FIG. 14: (a) The comparison between V(A, a, P, x) (solid 
line) and V^a;) (dashing line); (b)The propagation of the 
light beam with an input intensity profile described by 
\ipi{A, a,P,x)\ 2 . Here w = 10, fj, = 1/^3, A = 48.257, a = 
—0.348, P = 0.0141 and the degree of nonlocality Wo/(V3fi) = 
10. 



nonlocal case. Numerical simulations confirm that the 
soliton solutions in the 2nd perturbation can describe 
the fundamental and second order soliton states of the 
NNLSE in high quality. For the nonlocal case of the 
exponential-decay type nonlocal response, the Gaussian- 
function-like soliton solution can't describe the funda- 
mental soliton state of the NNLSE JJ) exactly even in 
the strongly nonlocal case, that greatly differs from the 
nonlocal case of the Gaussian function type nonlocal re- 
sponse. The functional dependence of the nonlocal soli- 
ton's power and phase constant on its beam width are 
greatly different from that of the local soliton. In the 
strongly nonlocal case, the soltion's power and phase con- 
stant are both in inverse proportion to the 4th power 
of its beam width for the nonlocal case of the Gaussian 
function type nonlocal response, and are both in inverse 
proportion to the 3th power of its beam width for the 
nonlocal case of the exponential-decay type nonlocal re- 
sponse. 



APPENDIX A: HOW TO CALCULATE THE 
PARAMETERS A, a, P WITH EQUATIONS. pOjl 

In principle, V(A, a, [3, x) and the parameters A, a, (3 
can be found by solving Eq. i(T§|> and Eqs. lj2"Ujl di- 
rectly, but these tasks are considerably involved. Here we 
present a fixed-point method to calculate these parame- 
ters A, a, (3 for a fixed value of fi. Firstly corresponding 
to V(A, a, 0, x) in Eq. iflUjl . we define 

U{a,/3,x) = — 2 . (Al) 

For an arbitrary pair of initial values of ckojA) with 
suitable order of the magnitude, we can calculate 
U(ao, (3q, x). Let 

Ai = \J M 4 (7< 2 )(«o,^5y' ( A2 ) 

ai = A?lrt 4 >(a ,/3b,0)/4!, (A3) 
(3 1 =A 2 1 U^(a (h 0o,O)/&.. (A4) 
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For such a pair of values of a%,Pi, we can find an- 
other U(ai , (3i , x) . Again we obtain another set of values 
{A 2 ,ct2, @2}- Repeating these steps of calculations, we 
can obtain series sets of values {A 2 , 0:3, P3}, {A3, a\, (3^}, 
and so on. The difference between {A m ,a m , /3 m } and 
{A m+ i, a m +i, Pm+i} will approaches zero as the num- 
ber of m approaches infinity. To some accuracy, we can 
calculate parameters A, a, (3 for a fixed value of /i. 
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For such a pair of values of we can find another 

[/(&). Again we obtain another set of values {A 2 , a 2 ,(3 2 }. 
Repeating these steps of calculations, to some accuracy 
we can calculate parameters A, a, (3 for a fixed value of 




FIG. 15: (a) The comparison between R(x) and 
\ip\(A, a, (3, x)\ 2 . Dashing line: R(x)/R(0); solid line: 
|^it4,a, A^lVlVMAa, AO. 65)j 2 ; (b)The propagation of 
the light beam with an input intensity profile described by 
\ipi(A,a,f3,x)\ 2 . Here w = 0.5, fj, = l/Vi,A = 6.303,« = 
-0.609, (3 = 0.0270 and the degree of nonlocality w /(V3[i) = 
0.5. 




APPENDIX B: HOW TO CALCULATE THE 
PARAMETERS A, a, (3 WITH EQUATIONS. (15511 

For a fixed value of /x and one suitable point xo ^ (in 
this paper we set Xq — 2), corresponding to V(A, a, j3, x) 
in Eq. (|19|l we define 



U{x) 



V(A,a,/3,x) 
A2 ' 



(Bl) 



For an arbitrary pair of initial values of ao, /3q with suit- 
able order of the magnitude, we can calculate U(x). Let 



A 1 = 



24[C/(a:o)-C/(0)]-9(7' 



'(x )x + U"(x )x^ ( B2 ) 



i2 7U' (x a )x„~-12[U(x„)-U(0)]-U" (x )xl /g^s 



i2 !7"(a;o)a;g+8[i7(a;o)-{/(0)]-5{/'(a;o)a;o ^g 4 ^ 




FIG. 16: (a) The comparison between V(A, a, /3, x) (solid 
line) and V(x) (dashing line); (b)The propagation of the 
light beam with an input intensity profile described by 
\ip 2 (A,a,/3,x)\ 2 . Here w = 2, fj, = l/y/E,A = 19.694, a = 
— 1.327, f3 — 0.0606 and the degree of nonlocality u>o/(VEfi) = 
2. 
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